% visualize.m
global xx yy
% 定义图形保存目录
dir1 = ['./figure/']; 

% 检查是否存在目录，如果不存在则创建
if ~exist(dir1, 'dir')
    mkdir(dir1);
end
nts=800
% 循环遍历所有时间步
for i = 1:nts
    if restart == 1
        i = nts_start + i;
    end
    
    % 加载当前时间步的数据
    load([saldir 'dat', sprintf('%4.4d', i), '.mat']);
    
    % 创建一个新的图形窗口
    figure('Name', sprintf('Time Step %d', i));
    set(gcf, 'Position', [500, 500, 1000, 500]); % 设置图形窗口的位置和大小
    
    % 绘制并保存温度分布图
    subplot(1, 2, 1);
    colormap jet;
    pcolor(yy(2:end-1), xx(1:end), squeeze(T0i(1:end, 2:end-1)));
    colorbar;
    shading interp;
    title(sprintf('T0 at t=%d', i));
    xlabel('y');
    ylabel('x');
    
    subplot(1, 2, 2);
    colormap jet;
    pcolor(yy(2:end-1), xx(1:end), squeeze(Ti(1:end, 2:end-1)));
    colorbar;
    shading interp;
    title(sprintf('Ti at t=%d', i));
    xlabel('y');
    ylabel('x');
    
    % 保存图形到文件，不显示图形窗口
    print(gcf, '-dpng', [dir1, sprintf('T0_Ti_%d.png', i)]);
    close(gcf); % 关闭图形窗口，避免闪烁
end